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Abstract 

Facilitated spin models on random graphs provide an ideal microscopic realization of the mode- 
coupling theory of supercooled liquids: they undergo a purely dynamic glass transition with no 
thermodynamic singularity. In this paper we study the fluctuations of dynamical heterogeneity 
and their finite-size scaling properties in the f3 relaxation regime of such microscopic spin models. 
We compare the critical fluctuations behavior for two distinct measures of correlations with the 
results of a recently proposed field theoretical description based on quasi-equilibrium ideas. We 
find that the theoretical predictions perfectly fit the numerical simulation data once the relevant 
order parameter is identified with the persistence function of the spins. 
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I. INTRODUCTION 



Facilitated spin models (FSM) have been proposed as an elegant way to rationalize the 
idea that the glass transition is purely dynamic in nature This contrasts the view 
that glassy phenomena can be understood in terms of a rugged free-energy landscape, as 
it happens in systems with complicated static interactions. Recent results have shown, 
however, that a common mechanism could be at the basis of the slowing down of relaxation 
both in FSM and in landscape dominated systems. Important insight has come from the 
study of cooperative FSM on random graphs, as they provide a suitable mean field limit in 
which the possibly sharp dynamic crossover observed in finite-size systems becomes an actual 



transition in the thermodynamic limit. In 



3| it has been shown that relaxational dynamics 



of FSM on random graph verifies the universal relations predicted by Mode Coupling Theory 



(MCT) of supercooled liquid 



iquids [4 



5[ and, more generally, that the formal structure of MCT 



singularities is closely related to that appearing in the bootstrap percolation analysis of 



FSM 



In some cases it is known that rugged landscapes and kinetic constraints provide 



a. dual description of slow relaxation [3, l£] The* recently studied ease of the XOR-SAT 



optimization problem of information theory {9] explicitly shows the emergence of a random 
first order transition scenario and the highly nonlocal nature of the relation between such 
dual descriptions. 

In this paper we address the question of whether the equivalence of description extends 
to dynamical heterogeneities and fluctuations. We focus, in particular, on the so-called (3 
relaxation regime, i.e. the time window in which the correlation function approaches its 
plateau value corresponding to the glass transition, and the system explores the configura- 
tion space near the ideal arrested state. We would like to test the predictions of a recent 
theory that describes fluctuations in the ft regime close to the MCT transition in terms of a 



cubic field theory in a random field (lOj. This theory, initially formulated for systems with 
quenched disorder, identifies the overlap with the initial condition as the order paramenter 
in terms of which a field theoretical description of fluctuations is possible. A priori, we do 
not expect such a theory to apply in systems in which the Hamiltonian is trivial (as in FSM) 
and the overlap does not give a sufficient determination of metastable states. 

The scheme of the paper is the following. In the second section we briefly review the 



theory of 



101 ] . In the third section we introduce the model and compare the numerical 



simulation results to the theory predictions for two different measures of correlations. In the 
last section we discuss the main implications of our findings and conclude the paper. 



II. FIELD-THEORETICAL APPROACH TO FLUCTUATIONS IN THE (3 REGIME: 
A MINI-REVIEW 



The theory proposed in 10| deals with equilibrium fluctuations of the local overlap, 
Qx{ s , s o), between an initial configuration, s = s(0), and the one visited by the system at 
a certain time t during its evolution, s = s(t). The appropriate definition of overlap may 
be system dependent and some coarse graining is generally assumed. For example, in liquid 
systems one can divide the space in small cells centered around the position x and compare 
the local particle density between the current and the initial configuration. In spin systems 
the most natural definition is the usual scalar product between the two spin configurations. 

The theory suggests to eliminate the temporal dependence of the relevant physical quan- 
tities from the description, starting from the observation that due to time scale separation, 
when the dynamical MCT transition is approached, dynamics can be described as a random 
walk between metastable states, equilibration within a state establishes before a new state 
is found. In such conditions one can eliminate time in favour of the local overlap with the 
initial state and describe fluctuations in a restricted equilibrium ensemble with probability 
measure: 



^-j exp y-PH{s) + Jdxv(x)q x (s,s Q )J , (1) 



where H(s) is the system Hamiltonian and = l/k B T is the inverse temperature. The 
Lagrange multiplier u(x) is fixed in such a way that the average values of the local overlap 
Qx{s, So) take some prescribed values close to the plateau value, q EA . The crucial assumption 
in the definition of the restricted measure, Eq. (pQ), is that the overlap provides a sufficient 
determination of the metastable states. 

A deep analysis of such a theory close to the ideal MCT transition (when activated 
processes are neglected) shows the emergence of a universal description in terms of the local 
overlap fluctuation 

<t>( x ) = Qx(s, s ) - g EA , (2) 
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which is encoded in the effective action: 

'1 



S 



dx 



(V0) + e<j)(x) + g<j>(x) d + h(x)(j)(x) 



(3) 



where e is a parameter quantifying the distance from the glass state (e.g. e = T — Ta) and 
g is a coupling constant. The "external field" h(x) account for the effects of heterogeneity 
in the initial condition: it is a random Gaussian variable with zero mean and variance 



E(h{x)h{y)) = A8{x-y). 



(4) 



The effective action can be used to study the finite-size corrections in mean-field glassy 
systems. In that case one only deals with global overlap fluctuations and so one just needs 
to replace the space integral with an overall volume factor N. The field variance A scales 
then as 1/N and the Gaussian theory suffices to understand finite-size scaling properties. 

In order to characterize fluctuations two types of four point correlation functions can be 
introduced: 



Xb*(t) = N (<Kt)> 2 - (<Kt)> 



(5) 
(6) 



We denote by the angular brackets, (■••), the average over trajectories that start from 
the same initial condition. This was called iso-configurational average in ll[ [l^]. In our 
stochastic dynamics the iso-configurational average represents the average over the thermal 
noise along the trajectories. The initial condition, denoted by s(0) = sq, will always be 
chosen as an equilibrium configuration and the corresponding average will be denoted by 
the overbar 7TT . Thus, Xth{t) characterizes thermal fluctuations for fixed initial condition, 
while Xhet{t) characterizes the fluctuations with respect to the initial condition. As typical in 
random field systems Xth an d Xhet have different scaling properties. If the quasi-equilibrium 
assumption of the theory is verified one can then eliminate the explicit time depedence of 
fluctuations in favor of 6, and it turns out that 



X th (4>,e,N) = m F th (<f>m,eN*), 
Xhet(0, e, N) = J- het (0iV^, eiV^) . 



(7) 
(8) 



The questions that we address in this paper are: i) does this description apply to FSM? ii) 
what is the relevant measure of correlations that characterizes metastability in FSM? 



Some of the basic ingredients of the theory, namely time scale separation and quasiergodic 
exploration of metastable states are manifestly true in FSM. What is not obvious in FSM 
is that the conventional spin overlap does provide a sufficient determination of metastable 
states. The analysis of FSM shows that the glass singularity is associated to a core or 
bootstrap percolation transition of frozen spins. The order parameter is therefore the fraction 
of frozen spins in the core. This suggests that the persistence function, p(t), i.e. the 
fraction of spins that have not flipped over the time interval [0, t], is the natural correlation 
function that provides the correct determination of metastable states. On the contrary, the 
conventional spin overlap, q(t), can be problematic in this respect, see e.g. the discussion in 
Ref. Ideally, we believe that it should be possible to write a MCT equation for p(t). Even 
more optimistically one could hope to write a field-theory like Eq. Unfortunately, this is 
a nontrivial task because p(t) is a highly non time-local quantity and its theoretical analysis 
is difficult. In this paper we limit ourself to discuss how the predictions, Eqs. (IT51) and (Tbfj) . 
are met if one concentrates on the persistence function and the conventional overlap. 



III. NUMERICAL SIMULATIONS OF FACILITATED SPIN MODELS ON RAN- 
DOM GRAPHS 

We have studied numerically the critical dynamics fluctuations of a facilitated system of 
N Ising spins, = ±1 with i = 1, . . . , N, on a regular random graph with coordination 
number z. The system has a trivial Hamiltonian 

h N 

i=l 

while spins evolve according to a Metropolis-like constrained dynamics: at each time step a 
randomly chosen spin, s«, is flipped with transition probability: 

w(si -»■ -Si) = min{l,e- hS!/fcBT } , (10) 

if and only if at least / of its z neighboring spins are in the state —1. Without loss of 
generality we set h = k B hereafter. We consider the cooperative case / = 2 with z = 4, 
which is characterised by a discontinuous glass transition at a critical value Ta ~ 0.48 
where the long-time limit of the persistence function, p(t), jumps from zero to a finite 
value p EA ~ 0.67. In full analogy with the mixed discontinuous-continuous nature of the 
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FIG. 1. Parametric plot of the two components of persistence fluctuations Xhet(t) an d Xth(t) vs - 
the average persistence, p(t), at temperature T = 0.5 close to the glass transition (T<j ~ 0.483) and 
for various values of the system size N. 



non-ergodicity parameter in MCT, the persistence has a square root singularity coming 
from low temperatures and share many other salient features of MCT 2], y], If], 13]. Such 
features are intimately accompanied by a substantial growth of dynamical heterogeneity as 
the glass transition is approached. Previous studies 14| have stressed the importance of 
the emerging time- reparametrization invariance of soft modes fluctuations during the aging 
dynamics 15l-ll7l|. Other manifestations of dynamic heterogeneity that have been studied 



include the Kovacs memory effect 



181 ] and the four-point dynamical susceptibility 



13 



19]. 



We shall focus here on the system size dependence of fluctuations at criticality. 

Apart from their conceptual underpinning such type of models offer the nontrivial advan- 
tage of requiring no thermal equilibration therefore allowing the investigation of relatively 
large system sizes and samples with a moderate computational effort. Efficient simulations 
are performed by means of a faster-than-the-clock algorithms 2^|. We have simulated spin 
systems of size ranging from iV = 128 to iV = 8192 for a number of samples of the order 
about 10 3 to 10 4 , depending on the system size. To measure Xhet and Xth we have generated, 
for each sample i — 1, . . . , M, two independent dynamical trajectories k — 1, 2 starting from 
the same equilibrium configuration and evolving with different thermal noise. Then for each 
sample and each trajectory we measure the persistence function 0^^, where the initial state 
is equal for the two dynamical trajectories. The two components of persistence fluctuations 
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FIG. 2. Scaling of persistence fluctuations in the (3 relaxation regime at the critical glass temper- 
ature T d = 0.483. 

are thus estimated as 




Notice that since quenched disorder is absent in such FSM there is no additional source of 
fluctuations apart from that one due to the topology fluctuations in the random generation 
of graphs. We have checked that such an effect is substantially negligible as expected from 
the self-averaging property of random graphs. 

In fig. 1 we illustrate how the two components of persistence fluctuations look like in the 
quasi-equilibrium description when they are expressed in terms of the persistence "clock" . 
We observe that heterogeneity fluctuations are almost one order of magnitude larger that 
thermal ones though they present qualitatively the same behaviour. Results are similar to 
those obtained in the p-spin spherical model. We then study the finite-size scaling of Xth 
and Xhet, which at the glass transition temperature, T = Td (e = 0), reads 

XM,0,N) = N$ T th (<f>N l *,0), (13) 
Xhet(0, 0, N) = N$ T het ((f)N*,0) , (14) 

where = p — p E A- We have first checked that the curves of XthN* and XhetN^ as a function 
of p obtained for various values of N cross precisely at the exactly known value of p EA . 
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FIG. 3. Heterogeneity and thermal overlap fluctuations at criticality T = (g EA ~ 0.4). 

In fig. 2 we see that, in agreement with the analytical predictions, the different curves of 
rescaled fluctuations plotted against the scaling variable N^cj) produces a pretty nice data 
collapse in a relatively wide region around the origin. We can thus summarise our first 
test by stating that the theory predictions are excellently verified when the correlations are 
measured in terms of the persistence function. 

Next, we study equilibrium fluctuations of the conventional overlap, q(t), that for our 
system should be more properly defined as 

"l - 

-2>(0) Si (f)-m(0)m(f) 



m 2 (0) 



i=i 



(15) 



where m is the magnetization density. Notice that q(t) exhibits features similar to the 
persistence function, i.e., for large system sizes its long-time limit jumps at from zero 
to q EA ~ 0.4 and has a square-root singularity below T d . It can play, consequently, the 
legitimate role of order parameter. However, this does not guarantee that it provides a 
sufficient determination of metastable states. Indeed, we find that two components of the 
overlap fluctuations do not obey the predicted finite size scaling, Eqs. (114"]) and (fT3|) . On 
one hand, we observe that the size dependence of heterogeneity fluctuations is substantially 
weaker than the expected a scaling form is recovered only if Xhet is rescaled by a factor 
ra t ner than N~^, see fig. 3a. On the other hand, we find that thermal fluctuations 
are of order one in the whole range q > q EA regardless the system size as it is observed in 
fig. 3b. These results confirm that although the overlap can be used as an order parameter 
to describe the glass transition, it yields only a partial characterization of metastability in 
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FSM [9]. Therefore, when the nature of metastable states in the system under consideration 
is not exactly known it is a rather delicate task to give a universal description of finite size 
critical fluctuations. 



IV. DISCUSSION AND CONCLUSION 

We have investigated the finite-size scaling properties of critical fluctuations in the /3 
relaxation regime of a cooperative facilitated spin model on a regular random graph (or Bethe 
lattice) for two different correlation measures. We tested our numerical simulation results 
against the predictions of a recent field theoretical approach that describes fluctuations in 
the (3 regime in terms of a cubic field theory in a random field. We find it remarkable 
that although FSM have been purposely devised to lack any thermodynamical content, the 
predictions of the quasi-equilibrium approach are exactly matched when correlations are 
measured in terms of the persistence function. In contrast, when correlations are measured 
through the conventional overlap function we find that theory predictions are not verified. 
Such a discrepancy can be traced back to the fact that the conventional overlap function 
does not provide a sufficient characterization of metastable states in FSM. Evidently, for 
finite dimensional systems, where the very notion of metastability is not sharply defined 
and requires the introduction of a suitable time scale possibly dependent on the observation 
time, the appropriate choice of relevant correlations becomes crucial. In the presence of 
an incomplete knowledge about the very nature of the metastability induced by the glass 
transition, some cautions is required when simulations or experiments are compared with the 
theory. For these reasons, it would certainly be interesting to extend the present investigation 
to finite dimensional FSM. This would give the possibility to assess, on one hand, the role of 
finite-size corrections to scaling, which in bootstrap percolation related problems are known 
to be anomalously large and, on the other hand, the relevance of the above scenario in 
systems in which the glass transition is avoided in the infinite size limit. Also, for continuous 
glass transitions, as those appearing in the simplest higher-order singularity scenario of MCT, 
we expect the above field theoretical description be still valid provided that the action S 
includes a quartic term, and the scaling behaviour of fluctuations is modified accordingly. 
Finally, from a more speculative theoretical perspective, we observe that since the aging 
dynamics of glassy systems can be interpreted in terms of quasi-equilibrium concepts one 
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would be tempted to generalise the field theoretical scenario we tested above in a suitable 
defined out of equilibrium regime of lenght and time scales. 
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